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Abstract 

Materials that combine magnetic and ferroelectric properties have generated in- 
creasing interest over the last few years, due to both their diverse properties and 
their potential utility in new types of magnetoelectric device applications. In this 
review we discuss recent progress in the study of such magnetoelectric multifer- 
roics which has been achieved using computational first-principles methods based 
on density functional theory. In particular, we show how first-principles methods 
have been successfully used to explain various properties of multiferroic materials 
and to predict novel effects and new systems that exhibit multiferroic properties. 
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1 Introduction 



Magnetoelectric multiferroics are materials which exhibit both magnetic order 
and ferroelectricity in the same phase. Such materials, although rare, have 
been known since the 1960s and various early review articles systematically 
classified their properties and behavior [1,2,3,4]. In the last few years the field 
of multiferroic materials has seen a tremendous boom, initiated in part by 
computational first-principles studies explaining the basic physics underlying 
their scarcity [5]. The subsequent flurry of research activity led to the discovery 
of several new multiferroic materials, as well as the identification of various 
strong coupling effects between their magnetic and ferroelectric degrees of 
freedom [6,7,8,9,10,11]. The prospect of using these coupling effects for new 
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types of devices in microelectronics and digital memory applications has in 
turn spawned further interest. 

This newer work on magnetoelectric multiferroics was recently reviewed in 
Ref. [12], with an emphasis on experimental and technological studies. The 
purpose of the present article is to provide an update on the recent activities 
in first-principles studies of magnetoelectric multiferroics, with the goal of 
illustrating the utility of these methods both in explaining unusual behavior 
in existing multiferroics, and in predicting new materials and novel effects. 
An extensive review of density functional theory-based studies of multiferroic 
materials before 2002 was given in Ref. [13]. Since then, first-principles studies 
of multiferroic materials have made considerable progress; here we summarize 
these recent developments and point out possible directions for future work. 

Ref. [13] reviewed two major questions which were successfully answered by 
first-principles electronic structure calculations. The first question, "Why are 
there so few magnetoelectric multiferroics?" [5] is now understood to result 
from a chemical incompatibility between magnetism and conventional ferro- 
electricity. This incompatibility is related to the fact that the most common 
mechanism for ferroelectricity in perovskite-structure oxides involves the pres- 
ence of a transition metal cation on the perovskite B site with a formal d° 
electron configuration. In fact, first-principles electronic structure calculations 
have played a major role in elucidating the role of the d° electron configu- 
ration in the ferroelectric instability in perovskite oxides [14,15,16]. On the 
other hand, for magnetism to occur in such transition metal oxides a partially 
filled d shell is indispensable. The existence of at least some known materials 
where magnetism and ferroelectricity coexist then poses the second question 
addressed in Ref. [13]: "Why are there any magnetoelectric multiferroics?" 
[17]. Clearly, if the ferroelectricity in a multiferroic is caused by a d° cation, 
a different cation is needed to introduce the magnetism; alternatively, if only 
magnetic cations with partially filled d shells are present, a new mechanism for 
ferroelectricity is required. To our knowledge, all known multiferroics adopt 
the latter scenario. An alternative mechanism for ferroelectric "off-centering" , 
which was discussed in Ref. [13], is the stereochemically active lone pair in 
cations such as Pb 2+ or Bi 3+ . Again, first-principles calculations were invalu- 
able in showing that this mechanism is responsible for the ferroelectricity in 
multiferroics such as BiMn03 or BiFe03 [18,19]. In these systems, the Bi 3+ 
cation has a formal valence electron configuration of 6s 2 6p° and the energy 
of the system can be lowered by off-centering of the Bi ion with respect to 
its oxygen surrounding, which leads to a hybridization of both Bi 6s and 6p 
states with O 2p orbitals and a localization of the lone pair on one side of the 
Bi ion [18]. 

The organization of this article is as follows. In the following section we give 
a brief overview of some methodological aspects of first-principles electronic 
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structure calculations that are especially relevant to the study of magnetoelec- 
tric multiferroics. Some recent methodological developments have contributed 
to the current boom in magnetoelectric multiferroics, since they facilitate a 
realistic description of the electronic structure of multiferroic materials and 
allow the calculation of quantities that are difficult to access experimentally. 

Following this brief methodological part, we discuss how first-principles calcu- 
lation have led to substantial progress in the understanding of multiferroic ma- 
terials during the last few years. We begin with cases in which first-principles 
calculations have explained properties of existing multiferroics. As examples 
we discuss the identification of the mechanism responsible for ferroelectricity 
in the hexagonal manganite YMnC>3, a material which is representative for a 
whole class of multiferroics, as well as the resolution of controversy surround- 
ing confusing and contradictory experimental data in multiferroic BiFe0 3 . 
Following that, we show how first-principles calculations have made impor- 
tant predictions in advance of their experimental observation. Here, we also 
discuss two examples: The identification of a mechanism for coupling between 
magnetic and ferroelectric order parameters that could facilitate the intrigu- 
ing effect of electric-field induced magnetization switching, and the ab initio 
design of new multiferroic materials with superior materials properties. Fi- 
nally, we give an outlook on future research in the field, and point out where 
first-principles methods can further contribute to the study of magnetoelectric 
multiferroics. 



2 Computational tools 

In this section we briefly summarize recent developments in electronic struc- 
ture methods which are particularly important for, and in one case were mo- 
tivated by, the study of multiferroic systems. For a complete description of 
density-functional-theory-based first-principles methods we direct the reader 
to one of the excellent texts which are available, e.g. Refs. [20,21]. 

The first electronic structure calculations of multiferroic systems were per- 
formed using the local spin density approximation (LSDA) to density func- 
tional theory [20]. Although the LSDA is a well-established technique that con- 
tinues to make invaluable contributions to our understanding of ferroelectrics 
(as well as many other materials), it has a number of limitations, among them 
its well-known underestimation of the size of the band gap for most insulat- 
ing materials. For systems containing localized electrons with strong Coulomb 
correlations such as transition metal oxides the use of the LSDA can be to- 
tally inadequate, since it can lead to metallic solutions for systems that are 
known experimentally to be insulators [22]. In the case of magnetic ferro- 
electrics this is particularly problematic, since it prohibits the calculation of 
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the spontaneous ferroelectric polarization and other indicators for structural 
instabilities, such as Born effective charges [23]. Therefore, the most important 
recent progress in first-principles methodology, with respect to the study of 
multiferroic materials, has been the development of methods that facilitate a 
realistic description of the electronic structure of magnetic systems containing 
strongly localized d or / electrons. 

Two methods that attempt to cure some of the deficiencies of the LSDA in 
the treatment of localized electrons were developed recently and have been 
applied to study multiferroic materials: the so-called "LSDA plus Hubbard 
U" method (LSDA+^7) [24] and a self-interaction corrected pseudopotential 
method ("pseudo-SIC") [25]. Although the physical motivation behind these 
two methods is different, they lead to very similar results in practice. 

The development of the pseudo-SIC method was actually motivated by the 
inability of the LSDA to produce a band gap for multiferroic YMJ1O3 (see 
section 3.1 below). In general, self-interaction corrected (SIC) methods are 
based on the well-known fact that the failures of the LSDA can be, at least in 
part, attributed to the presence of the self-interaction (SI) in the LSDA energy 
functional, that is the interaction of an electron charge with the Coulomb and 
exchange-correlation potential generated by the same electron [20]. The SI is 
strong for spatially localized electrons such as 3d and 4/ electrons, and leads 
to an underestimation of the binding energies, on-site Coulomb energies and 
exchange splittings of the corresponding states, while the hybridizations of 
cation d and anion p states and the corresponding band widths are usually 
overestimated. Possible strategies for eliminating the SI in density functional 
theory are long-standing issues, but often require a large computing effort 
even for materials with small unit cells. In the pseudo-SIC method [25], the SI 
correction is calculated for the corresponding atoms (as in Ref. [26]), and is 
then scaled by the electron occupation numbers of specific atomic-like orbitals 
calculated self-consistently within the crystal environment. This approach can 
be efficiently implemented within a standard pseudopotential method and 
allows the SI originating from localized, hybridized, or completely itinerant 
electrons to be discriminated, which permits the treatment of insulating as 
well as metallic compounds, with minimal computational overhead beyond 
the LSDA [25]. 

The motivation behind the development of the LSDA+[7 method is some- 
what different, although in practice this method leads to very similar results 
as the pseudo-SIC method. In the LSDA+[7 method the Coulomb interaction 
between specific localized orbitals is explicitly taken into account by adding 
a term to the total energy that resembles the well-known "Hubbard Hamil- 
tonian" used to describe Mott-Hubbard metal-insulator transitions [27]. The 
resulting potential leads to an energy splitting between occupied and unoccu- 
pied states. The Hubbard parameter, U, which describes the strength of the 
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screened electron-electron interaction, can in principle be calculated ab initio 
but in practice it is often treated as an empirical parameter. Simultaneously, 
a "double-counting correction" is subtracted from the total energy in order to 
account for those interaction effects between the localized electrons that are 
already included in the conventional LSDA. Different forms for this correction 
term are in use, some of which are in part also designed to make the resulting 
energy functional exhibit certain desirable properties [28]. These ambiguities 
and the partly empirical construction of the Hubbard correction is somewhat 
unsatisfying and the resulting method is strictly speaking not purely first prin- 
ciples. Nevertheless, the LSDA+[7 method is nowadays used more and more 
routinely in electronic structure calculations for d and / electron systems and 
has undoubtedly led to considerable progress in the field. 

Simultaneously with the tremendous improvements in the first-principles de- 
scription of magnetic insulators provided by the LSDA+U and pseudo-SIC 
methods, the so-called "modern theory of polarization" has revolutionized 
the first-principles study of ferroelectric materials [29,30,31]. This theory pro- 
vides both conceptual understanding of the spontaneous polarization in bulk 
crystals, and computational tools to calculate it from the Kohn-Sham wave- 
functions. Complications in calculating the electric polarization in bulk crys- 
tals arise from the periodic boundary conditions employed in most solid state 
calculations. This periodicity leads to an ambiguity in the value of the polar- 
ization modulo an integer multiple of the so called "polarization quantum", 
such that 



where P and P' are two allowed values of the polarization, n is an integer 
number, e is the electronic charge, V is the unit cell volume, and R is a 
primitive lattice vector of the underlying Bravais lattice. The polarization of 
an infinite periodic crystal is therefore represented by a lattice of values, the so- 
called "polarization lattice" . In an experiment only polarization changes can 
be measured and the corresponding value has to be calculated as a difference 
between two points of the polarization lattices corresponding to the initial 
and final states of the system. The corresponding change in polarization is 
uniquely defined as long as the system stays insulating for all intermediate 
states transforming the initial to the final state. 

Figure 1 shows the evolution of the polarization lattice in multiferroic BiFeC>3 
calculated for different points along a ferroelectric switching path leading from 
the negatively polarized ferroelectric state ( — 100% distortion), through the 
corresponding centrosymmetric structure (0% distortion), to the positively 
polarized state (+100% distortion) [19]. The vertical columns of black circles 
represent various points of the polarization lattice for a fixed amount of dis- 
tortion, separated by the polarization quantum of 185.6 /iC/cm 2 . The dashed 
lines show how the polarization lattice changes along the path. 
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Fig. 1. Change in polarization P along a ferroelectric switching path for multiferroic 
BiFeOs. The possible values of P for fixed distortion differ by multiples of the 
polarization quantum, here 185.6 //C/cm 2 . From Ref. [19]. Copyright (2005) by the 
American Physical Society 



The polarization lattice can be calculated as a gauge invariant geometric phase 
of the wave-functions, the so-called "Berry-phase" [29,30,31]. This expression 
solves the problem arising from the unboundedness of the quantum mechanical 
position operator with the infinitely extended periodic wave-functions in the 
solid and is compatible with the "multivaluedness" of the polarization; the 
ambiguity of the polarization modulo a quantum translates into an ambiguity 
of the Berry-phase modulo integer multiples of 2n. 



In addition to "beyond LSDA" methods like pseudo-SIC and LSDA+L^ which 
result in stable insulating solutions for most magnetic oxides, and the modern 
theory of polarization, which allows an accurate calculation of experimentally 
measurable polarization differences, there are many other developments in 
first-principles-based methods that are relevant to the study of multiferroic 
materials, which we mention briefly. For example, many multiferroic systems 
have complex noncollinear spin structures, and coupling effects between the 
polar and magnetic degrees of freedom can be due to spin-orbit coupling. 
Therefore, a fully noncollinear treatment of the spinor wave-functions and the 
inclusion of spin-orbit coupling in the calculations is essential to describe many 
of the interesting properties of multiferroic materials. Also, the analysis of fer- 
roelectric instabilities using calculated phonon spectra is a very powerful tool 
to explain the appearance of multiferroic properties and predict the existence 
(or absence) of such instabilities in designing new magnetoelectric multiferroic 
materials. 
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3 First-principles explanation of experimental observations 

We now outline recent first-principles work that has led to substantial progress 
in the understanding of known multiferroic materials. First we describe stud- 
ies aimed at understanding the unusual mechanism for ferroelectricity in the 
hexagonal manganite materials. Then we discuss a range of seemingly contra- 
dictory experimental reports of both the magnetic and ferroelectric data in 
multiferroic BiFeC-3, and show how first-principles calculations were able to 
resolve some of the controversy. 



3.1 Elucidating the origin of ferroelectricity in hexagonal manganites 

Hexagonal YMn0 3 is the prototype of the class of isostructural multifer- 
roic materials with the chemical formula i?Mn03 (i?=Ho-Lu, Y) and room 
temperature space group PQ^cm. These materials are ferroelectric with high 
ferroelectric Curie temperatures of ~ 1000 K, and exhibit noncollinear an- 
tiferromagnetic order below T N ~ 100 K. Although there have been many 
recent experimental studies of the corresponding rare-earth systems, primar- 
ily HoMn03, which show a complicated phase diagram at low temperatures, 
the theoretical first-principles work has focused on YMnOs, partly because of 
the absence of 4/ electrons in this system, which are difficult to treat using 
pseudopotential-based methods. The main effort has been directed to under- 
standing the mechanism for the occurrence of ferroelectricity, which is clearly 
of an unconventional origin, since the Mn 3+ cation does not have a d° electron 
configuration and Y 3+ does not have a lone pair, so the system contains no 
"ferroelectrically active" ions. Initial progress in understanding the ferroelec- 
tric behavior was hampered by the severe underestimation of the gap within 
LSDA [16], and by the fact that early experiments incorrectly reported an off- 
centering of the Mn ions from the centers of their polyhedra in the ferroelectric 
phase [32]. 

The first important contribution from first-principles calculations was the 
analysis of the electronic structure of hexagonal YMnOa using the LMTO- 
ASA method [33]. These calculations, using the early experimentally reported 
paraelectric and ferroelectric structures, showed that the Mn 4+ ion in hexag- 
onal YMn03 is not a Jahn- Teller ion despite of its d A electronic configuration, 
because the trigonal crystal field splitting in this material does not lead to 
an orbital degeneracy [33]. The width of the band gap was shown to be very 
sensitive to the magnetic ordering, increasing (within the LSDA) from about 
eV for ferromagnetic [33] and A-type antiferro magnetic [34] order to 0.47 eV 
for a more complicated but still collinear antiferromagnetic spin arrangement 
[33]. The width of the gap was also shown to strongly increase if correlations 
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were incorporated in the calculations through an LSDA+£7 treatment [34,33]. 
This led to the conclusion that both magnetic ordering and correlation effects 
have a strong influence on the electronic structure, and play an important role 
for the gap formation in hexagonal YMn0 3 . These calculations also showed, 
however, that the overall electronic structure is not strongly affected by the 
experimentally reported structural distortions leading from the high temper- 
ature paraelectric phase to the low-temperature ferroelectric phase. The only 
significant change observed in the calculations was a slightly stronger p-d hy- 
bridization due to the off-center displacement of the Mn ion within its sur- 
rounding oxygen cage that was reported in the early experiments [32]. 

Subsequent plane-wave pseudopotential calculations confirmed that such dis- 
placements of the Mn cation along the hexagonal axis indeed lead to a certain 
degree of rehybridization between O 2p and the Mn 3d z 2 orbitals which are ori- 
ented along the polar c direction, whereas displacements within the c plane are 
not accompanied by large changes in chemical bonding [35,16]. These obser- 
vations, together with the experimental reports of Mn c-axis displacements, 
led to the suggestion that a "one-dimensional a?°-ness", i.e. the unoccupied 
nondegenerate d z i orbital oriented along the displacement direction, could be 
responsible for the stabilization of the ferroelectric state [16]. Nevertheless, an 
energy-lowering due to the displacement of the Mn ions was not found in the 
first-principles calculations. It was unclear at this stage whether the absence of 
energy lowering was an artifact of the use of the LSDA, which does not lead to 
an appropriate description of the electronic structure of YMnC>3 (in fact, in the 
plane-wave pseudopotential calculations the system remained metallic even in 
the case of A-type antiferromagnetic ordering), or an indicator that the early 
experiments reporting a Mn c-axis displacement were incorrect. Therefore, a 
conclusive picture of the ferroelectricity could not be achieved. 

The first breakthrough came through a collaboration that combined progress 
on both the experimental and the theoretical side. First, new single crystal 
x-ray diffraction experiments, using reflections of the entire Ewald sphere, did 
not find an off-centering of the Mn ions within their oxygen cages [36]. Instead, 
the experimental data indicated that the Mn ions remain in the center of their 
oxygen polyhedra which instead undergo a collective rotation that leads to unit 
cell tripling. In addition, the Y ions move along the hexagonal c axis, and fill 
the space that is freed due to the rotation of the oxygen bipyramids (see Fig- 
ure 2). Second, on the theoretical side the development of the self-interaction 
corrected pseudopotential (pseudo-SIC) method [25], motivated by the inad- 
equacies of the LSDA, enabled a full structural optimization of YMn0 3 while 
giving a realistic description of the underlying electronic structure. These cal- 
culations confirmed the new experimental data, giving structural parameters 
very close to those obtained experimentally. The first-principles calculations 
also showed that the structural distortions in the low temperature phase do 
not lead to any significant charge redistribution. This is reflected in the values 
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Fig. 2. High temperature centrosymmetric (left) and low-temperature ferroelectric 
structure (right) of YMn03. In the low-temperature structure the oxygen polyhedra 
(turquoise) undergo a collective rotation that leads to a unit cell tripling. This unit 
cell tripling is accompanied by displacements of the Y cations (blue) along the 
(0001) directions which lead to an electric dipole moment. 

of the Born effective charges, which were found to be very close to the formal 
charges. It became clear that rehybridization is not responsible for stabilizing 
the ferroelectric state in YMnC>3 and that in fact the interplay between polar 
and nonpolar structural modes is essential to the ferroelectricity in this system 
[37]. 

A concise picture of this interplay between polar and nonpolar structural 
modes was obtained very recently by calculating and analyzing the phonon 
instabilities in the high temperature centrosymmetric structure [38]. It was 
shown that the only unstable phonon mode is a unit-cell tripling mode, which 
results in the experimentally observed collective tilting of the oxygen polyhe- 
dra and a buckling of the Y (0001) planes (K 3 mode). Although this mode on 
its own does not lead to an electric dipole moment, it reduces the symmetry 
of the system to the polar space group P6 3 cm, which corresponds to the ex- 
perimentally observed low-temperature structure. If the distortions that lead 
from the high-temperature centrosymmetric to the low-temperature ferroelec- 
tric phase are decomposed into contributions from different phonon modes, the 
contribution of this K% mode makes up for more than 80 % of the final distor- 
tion. The second largest contribution (about 15 %) is due to a zone-center po- 
lar mode (rj mode) which is compatible with the symmetry reduction due to 
the unstable K 3 mode but is stable in the centrosymmetric high-temperature 
structure. By mapping out the energy surface around the prototype phase 
as a function of the K 3 and the polar mode it was shown that the polar 
mode remains stable for all values of K 3 but that due to a symmetry allowed 
coupling term between the and K 3 modes the equilibrium value of the T^" 
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distortion amplitude is shifted to a nonzero value. The K% mode therefore acts 
as a "geometric field" on the mode. This behavior is indicative of an im- 
proper ferroelectric where the polar mode is not the primary order parameter 
and this mechanism has been suggested to be the driving force behind the 
ferroelectricity in YMnC>3 [38]. 

The example of YMnO"3 shows how first-principles calculations of phonon spec- 
tra and analysis of the calculated electronic structure can give a good explana- 
tion of the mechanism driving the ferroelectric phase transition in multiferroic 
materials. It is also apparent that the use of the LSDA+[7 and pseudo-SIC 
"beyond LDA" methods, and the resulting realistic description of the elec- 
tronic structure, has been essential for such progress to be achieved. 



3.2 Resolving the true magnetization and polarization in BiFe0 3 

The determination of the spontaneous polarization and magnetization in BiFeC>3 
is an example of how valuable first-principles calculations are for determining 
intrinsic material properties and explaining contradictory experimental data. 
In a real sample, effects such as epitaxial strain, defects, microstructure, or 
interface properties can alter the physical properties of a material. Therefore it 
is sometimes extremely difficult to determine the intrinsic material properties 
unambiguously using only experimental data. First-principles calculations on 
the other hand can deliver reference data corresponding to a perfect sample 
containing absolutely no defects. After establishing the intrinsic value of a 
certain quantity, these calculations can also be used to determine how effects 
such as epitaxial strain, various kinds of defects, or the presence of interfaces 
and surfaces change this intrinsic value. 

Perovskite-structure BiFeC>3 is one rare example of a magnetoelectric multi- 
ferroic that exhibits both magnetic and ferroelectric ordering above room tem- 
perature. The bulk material becomes ferroelectric below Tc ~ 1103 K [39], 
adopting the R3c structure shown in Fig. 3. The magnetic moments of the Fe 
cations order antiferromagnetically (G-type) below T N 643 K [40,41], and in 
addition it exhibits a long-wavelength spiral-spin structure with a wavelength 
of about 620 A [42] . Until recently, the value of the polarization in BiFeO"3 was 
believed to be rather small. The only current-voltage measurements performed 
on BiFe0 3 single crystals reported a spontaneous polarization of 3.5 /zC/cm 2 
along the [100] direction (corresponding to 6.1 /iC/cm 2 along the polar [111] 
axis), but the corresponding hysteresis loops were far from being saturated 
[39]. Similar small values of the polarization were subsequently reported for 
(Bi . 7 Bao.3)(Fe . 7 Tio.3)0 3 films on Nb-doped SrTi0 3 (100) [43] and for poly- 
crystalline BiFeO"3 films [44] . However, all these measurements were hampered 
by the high conductivity of the available samples. 
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Fig. 3. BiFeC>3 crystallizes in a distorted perovskite structure with space group R3c. 
The cations are displaced from their ideal positions relative to the anions along the 
[111] direction. In addition, the oxygen octahedra are rotated around the [111] 
direction, alternately clockwise and counterclockwise (Bi: blue, Fe: green, O: red). 



Interest in BiFeOa has grown considerably over the last few years, following a 
report of large electric polarization and substantial thickness-dependent mag- 
netization in epitaxial thin films [7]. The large ferroelectric polarization and 
macroscopic magnetization were unexpected, since all earlier studies had re- 
ported only a small value of the polarization and no magnetization. Therefore 
the large polarization measured in the epitaxial thin films of BiFeOs was ini- 
tially ascribed to the effect of the heteroepitaxial constraint and the resulting 
change in the lattice parameters. Since a multiferroic with large polarization 
and magnetization above room temperature would be very desirable for tech- 
nological applications [12], a flurry of experimental activity ensued but could 
not resolve the confusion, with a range of reported polarizations ranging be- 
tween 2.2 yuC/cm 2 [44] and 158 ^C/cm 2 [45]. 

The spread in magnetization and polarization data therefore posed a number 
of questions: 

(1) What are the intrinsic values of the spontaneous polarization and the 
magnetization in BiFeOa? 

(2) How are the polarization and the magnetization in BiFe03 affected by 
epitaxial strain? 

(3) How are the polarization and magnetization affected by defects? 

In the following we show how all these questions have been answered using 
first-principles calculations. 
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3.2.1 Intrinsic polarization of BiFeOs 

In Ref. [7] the large polarization found in the epitaxial thin films of BiFe0 3 
was initially ascribed to the effect of the heteroepitaxial constraint and the 
resulting change in the lattice parameters. To test this hypothesis, the po- 
larization was calculated from first-principles, both in the bulk rhombohedral 
structure, and a hypothetical tetragonal structure, with its in-plane lattice pa- 
rameters constrained to that of the substrate. These first calculations seemed 
to corroborate the hypothesis that the large polarization was due to the epi- 
taxial constraint, since the calculated values for the two different structures 
were considerably different. However, the LSDA was used in these calculations 
and no insulating centrosymmetric reference structure (see Section 2) could 
be obtained. Therefore, only the absolute values of the polarizations were re- 
ported (that is the values at +100% and —100% distortion in Fig. 1) and not 
the differences in polarization between the ferroelectric structure and a cen- 
trosymmetric reference structure, which would correspond to the polarization 
measured in a ferroelectric switching experiment. Note also, that subsequent 
experiments showed that the structure in the films is in fact similar to the 
rhombohedral bulk structure, albeit with an additional monoclinic distortion 
in the case of a (001) oriented substrate [46,47,48]. 

Subsequent first-principles calculations using the LSDA+^7 method were able 
to obtain an insulating solution for BiFe0 3 in the centrosymmetric cubic per- 
ovskite structure and in the R3c structure, both of which are reasonable ref- 
erence structures for calculating the spontaneous polarization of BiFe0 3 in 
the rhombohedral bulk structure [19]. Calculation of the polarization change 
along the so-defined switching path resulted in a value for the spontaneous 
polarization in bulk BiFe0 3 of about 95 /uC/cm 2 along the [111] direction, 
consistent with the values measured in the thin films [7,46], but exceeding 
the measured values in the bulk samples by one order of magnitude. This re- 
sult suggested that the small values of the polarization reported in the past 
were in fact due to incomplete switching and the resulting unsaturated hys- 
teresis loops represented only a fraction of the true polarization of BiFe0 3 . 
The large polarization found in the thin films is then the first really saturated 
measurement of the full polarization of BiFe0 3 . 



3.2.2 Effects of strain and defects on the spontaneous polarization in BiFe0 3 

It is known from materials such as BaTi0 3 and SrTi0 3 that the epitaxial 
strain present in thin film samples can lead to drastic changes in the ferroelec- 
tric polarization compared to the corresponding bulk materials [49,50]. The 
good agreement between the calculated bulk value of the polarization and the 
experimentally observed thin film value in BiFe0 3 therefore seems to be in 
conflict with the strong strain dependence of polarization expected in general 
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for ferroelectrics. The lattice mismatch between BiFe0 3 and SrTi0 3 , the sub- 
strate material used in Ref. [7], amounts to about 2% compressive strain in 
the BiFe0 3 film, which is a substantial value. 

To address the epitaxial strain dependence of the ferroelectric polarization in 
rhombohedral BiFe0 3 a detailed first-principles study was carried out. Epi- 
taxial strain was imposed by constraining the in-plane lattice parameters over 
a range of non-equilibrium values, and allowing the structure to relax out-of- 
plane at each constraint [51]. In contrast to similar calculations for BaTi0 3 
and PbTi0 3 [52,53,54], it was found that the electric polarization in BiFe0 3 is 
not significantly affected by the presence of epitaxial strain. Even strain values 
as large as ± 3% lead only to changes in the polarization of a few percent, pro- 
vided that the material does not undergo a structural phase transition. (For 
comparison, the corresponding change in BaTi0 3 is more than 100% [52]). 
The relative strain-independence of the polarization in BiFe0 3 supports the 
notion that the large polarization found in the thin films is essentially equal 
to the bulk value and that the old measurements severely underestimated the 
polarization in BiFe0 3 . 

Recently, possible reasons for this different strain behavior between BiFe0 3 
and ferroelectrics like BaTi0 3 or PbTi0 3 have been addressed using first- 
principles methods [55]. It was shown that the variation of the spontaneous 
polarization with epitaxial strain in thin film ferroelectrics can be understood 
in terms of the piezoelectric and elastic constants of the corresponding bulk 
materials, and that a large strain dependence is not a common feature of all 
ferroelectrics. 

Ref. [51] also investigated the effect of oxygen vacancies on the spontaneous 
polarization in BiFe0 3 using first-principles techniques. It was shown that the 
incorporation of such vacancies leads to small relaxations around the vacancy 
site and to the formation of Fe 2+ next to the vacancy site, but that the value 
of the spontaneous polarization is not significantly affected by the presence of 
oxygen vacancies. 



3.2.3 Magnetic properties of BiFeO^ 

First-principles calculations have also been applied to understanding the mag- 
netic properties of BiFe0 3 , in particular the observation of magnetization in 
BiFe0 3 thin films and its absence in bulk samples. One possible source of mag- 
netization in a primarily antiferromagnetic material is weak ferromagnetism, 
i.e. a canting of the mainly antiferromagnetically oriented magnetic moments 
resulting in a small magnetization [56]. In bulk BiFe0 3 , the G-type antiferro- 
magnetic ordering [41] is modulated by a long wavelength spiral spin structure 
[42] which would lead to a cancellation of the magnetization on a macroscopic 
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scale. However, it was suggested that the thin film morphology might suppress 
the spiral spin structure so that the weak ferromagnetism becomes observable 
in thin films [56,57]. First-principles calculations indeed showed that a local 
canting does occur, and that its origin is the Dzyaloshinskii-Moriya interaction 
[58]. The magnitude of the calculated value for the magnetization is approxi- 
mately 0.1 lib/ (unit cell), in good agreement with recent measurements [57,59]. 
Therefore, this supports the interpretation proposed above for the origin of 
the magnetization in thin films of BiFeOa. 

The thin film measurements of Ref. [7] also reported a strong increase in 
magnetization with decreasing film thickness, reaching values of about 1 /xe/Fe 
for film thicknesses below 100 nm. Since neither calculations on ideal systems 
[56] nor recent measurements on highly stoichiometric BiFeC>3 films [59] were 
able to reproduce an increased magnetization, it was suggested that oxygen 
vacancies and the resulting presence of Fe 2+ could lead to the observed increase 
in the magnetization of very thin films [60]. However, supercell calculations 
corresponding to different concentrations of oxygen vacancies were not able to 
confirm this suggestion [51]. In these calculations it was found that removing 
a neutral oxygen atom from the supercell resulted in the formation of Fe 2+ 
on the sites adjacent to the vacancy site. The magnetic moments of such a 
pair of Fe 2+ cations are coupled antiferromagnetically and therefore do not 
result in a net magnetization. Also, the canting of the magnetic moments was 
not significantly affected by the presence of the oxygen vacancies, so that no 
increase in the macroscopic magnetization due to oxygen vacancies could be 
found. Ref. [51] also showed that the canting, and with this the magnetization, 
is not significantly enhanced by epitaxial strain. The reason for the increase 
in magnetization with decreasing film thickness in Ref. [7] is therefore still 
unclear, but from the first-principles calculations it can be concluded that the 
observed effect must be caused by effects other than epitaxial strain or neutral 
oxygen vacancies. 



3.3 Other systems 

Very recently, phonon energies were calculated for the spinel system CdCr 2 S4 
[61]. This material, which is known to be a rare example of a ferromagnetic in- 
sulator with a rather high magnetic Curie temperature, was recently reported 
to exhibit interesting coupling between dielectric and magnetic properties as 
well as a small electric polarization at low temperatures [62] . It was suggested 
that the relaxor-like behavior could be due to frustrated soft-mode behavior 
and that the ferroelectricity at low-temperatures could originate from a dis- 
placement of the Cr 3+ ion [62]. The first-principles calculation on the other 
hand showed that all phonon modes in this material are rather stable and that 
the scenario of a Cr 3+ ion moving off-center is extremely unlikely [61]. More 
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recent experimental results indicate that the observed behavior is probably 
defect-related and not an intrinsic effect since the observed features in the di- 
electric constant vanish after annealing the corresponding samples [63]. This 
represents another example how first-principles calculations can be used to 
separate intrinsic material properties from extrinsic effects like microstructure 
or defects. 

3.4 Summary of first-principles explanation of experimental observations 

It is clear from the above examples that first-principles calculations have made 
many contributions to our understanding of multiferroic phenomena, and that 
the interplay between theory and experiment has been invaluable for progress 
in the field. The example of YMnC>3 showed how the mechanism underlying 
unconventional ferroelectric phase transitions can be elucidated using first- 
principles techniques, and the example of BiFeC>3 showed how first-principles 
calculations can be used for resolving conflicting experimental results and to 
separate intrinsic materials properties from extrinsic effects such as defects or 
microstructure. 



4 First-principles predictions of new materials and novel phenom- 
ena 

In this section we change our focus compared to the previous section and dis- 
cuss examples where predictions have been made using first-principles tech- 
niques prior to their experimental verification. 

4-1 Electric-field induced magnetization switching 

The most important question arising from the presence of both spontaneous 
electric polarization and macroscopic magnetization in a multiferroic is whether 
these two properties are coupled, and, if so, in what manner. A coupling that 
facilitates the reversal of the magnetization by an electric field or conversely, 
the reversal of the electric polarization by a magnetic field, would be of par- 
ticular interest. In the following we show how first-principles calculations have 
been used to study this question, and to develop theoretical models for cou- 
pling between magnetic and structural/dielectric properties. 

Again, BiFeC>3 has been used as the prototype material for first-principles 
study. Calculations of the coupling between magnetism and structural distor- 
tions in BiFe0 3 showed that the Dzyaloshinskii-Moriya interaction [58], which 
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Fig. 4. Possible realization of electric-field induced magnetization switching in a 
weak ferromagnetic ferroelectric: if the polarization P is reversed, the canting of the 
antiferromagnetic sublattices, and with this the magnetization M, is also reversed. 



leads to the canting of the magnetic moments and the appearance of a macro- 
scopic magnetization, is indeed a result of the symmetry lowering structural 
distortions in this material (in the high symmetry cubic phase, weak ferromag- 
netism is not allowed by symmetry) [56] . In the case of BiFeOs there are two 
different symmetry lowering distortions that together reduce the symmetry 
from the cubic perovskite structure (Pm3m) to the ferroelectric rhombohe- 
dral ground state structure (R3c). One is a nonpolar R mode that leads to a 
unit cell doubling compared to the simple perovskite structure and rotations of 
the oxygen octahedra around the rhombohedral axis, and the other is a polar 
T mode that displaces all ions relative to each other along the same rhom- 
bohedral axis (see Figure 3). It was shown that the Dzyaloshinskii-Moriya 
interaction in BiFeOa is caused by the nonpolar mode and that reversal of 
that mode also inverts the magnetization direction [56]. 

This has important consequences for the possible realization of electric-field- 
induced magnetization switching, which was believed to be impossible due to 
symmetry reasons [64]. From the results obtained in Ref. [56] it became clear 
that such an effect is in principle possible in weak ferromagnetic systems in 
which the polar distortion, which results in ferroelectricity, also leads to the 
appearance of a Dzyaloshinskii-Moriya interaction, which results in weak fer- 
romagnetism. In such a scenario, the magnetization will be strongly coupled 
to the ferroelectric polarization and reversal of the ferroelectric polarization 
by an electric field will be accompanied by reversal of the macroscopic mag- 
netization [56] (see Figure 4). The search for a material where this effect can 
be realized is a topic of current research activity [65]. 
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4-2 Computational materials design 



The examples discussed in the previous sections illustrate that first-principles 
calculations are capable of correctly predicting numerous properties of mag- 
netoelectric multiferroics. This opens up the possibility of ab initio computa- 
tional design of entirely new multiferroic materials, with improved properties, 
in advance of their synthesis. Such a materials design approach was suggested 
in Ref. [66] and has been applied since then to search for new multiferroics 
with properties superior to those of existing materials. 

The scarcity of magnetoelectric multiferroics has rendered such an approach 
extremely desirable. For example, most magnetoelectric multiferroics are ei- 
ther antiferromagnets or exhibit only weak ferromagnetism, often below room 
temperature. On the other hand, for technological applications a multiferroic 
with large polarization and large magnetization, well above room temperature, 
and preferably with strong coupling between the two properties is needed. Cur- 
rently, no magnetoelectric material is known that fulfills all these requirements. 
Since the calculation of the properties of a variety of candidate materials is of- 
ten less costly than the corresponding experimental synthesis and subsequent 
characterization, a first-principles design approach is very appealing. 

Probably the first example for the ab initio design of a magnetoelectric mul- 
tiferroic was the prediction of ferroelectricity in BiMnC>3. This material was 
known to be ferromagnetic at low temperatures [67], and the calculation of un- 
stable zone-center phonon modes in cubic BiMnOs suggested the simultaneous 
existence of ferroelectricity [68]. Motivated by this prediction, experimental in- 
vestigations were performed and ferroelectric hysteresis loops were reported 
for polycrystalline samples as well as for thin films of BiMnC>3 [69] . The ferro- 
magnetism in BiMn0 3 is interesting in its own right, and is believed to be due 
to orbital ordering caused by a cooperative Jahn- Teller distortion [70]. Such an 
orbital ordering was reproduced recently by first-principles calculations using 
the experimentally obtained structure [71]. 

A first-principles search for multiferroism has also been performed for the 
closely related compound BiCr0 3 [72]. Prior to the theoretical study, little was 
known experimentally about the material; early papers reported antiferromag- 
netism with a weak ferromagnetic moment below 123K, and a structural tran- 
sition between pseudomonoclinic and pseudotriclinic perovskite structures at 
400K [73] . The first-principles calculations found an antiferromagnetic ground 
state, in agreement with the Goodenough/Kanamori/Anderson rules for su- 
perexchange interactions in magnetic insulators [74,75], and the calculation 
of the full phonon spectrum revealed several unstable phonon branches over 
the whole range of the Brillouin zone. The most unstable modes were found 
at the R and M points of the Brillouin zone corresponding to antiferroelec- 
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trie or antiferrodistortive phonon modes. Since the antiferromagnetic ordering 
in BiCrC>3 does not lead to a macroscopic magnetization, it was concluded 
that this material will not be of great technological relevance and a full struc- 
tural optimization of the ground state structure was not attempted [72]. A 
subsequent experimental study [76] confirmed the weak ferromagnetism, and 
observed a non-polar to polar structural phase transition, suggestive of ferro- 
electricity, at 440K. 

The difficulty of incorporating a substantial magnetization in a magnetoelec- 
tric multiferroic results from the dominance of magnetic superexchange in such 
transition metal oxides, which in most cases leads to antiferromagnetic near- 
est neighbor coupling. Although in certain cases the superexchange mechanism 
can lead to ferromagnetic coupling, the strength of such a ferromagnetic cou- 
pling is usually substantially weaker than in the antiferromagnetic case [75], 
resulting in low magnetic ordering temperatures. Therefore a recent proposal 
to exploit the strong antiferromagnetic superexchange and introduce macro- 
scopic magnetization by combining two magnetic ions with different magnetic 
moments, resulting in a ferrimagnetic spin-arrangement [77] is very promising. 
The ordered double perovskite system Bi 2 FeCr0 6 was suggested as a candi- 
date material, since the stereochemically active Bi 3+ ion should introduce a 
potentially ferroelectric structural distortion and the closed-subshell d 5 Fe 3+ 
and d 3 Cr 3+ ions should result in insulating behavior. 

Indeed, a ferroelectric ground state with space group R3 was found by struc- 
tural optimization using first-principles calculations within the LSDA+[7 method 
[77] . The corresponding structure is closely related to the R3c structure found 
in the well-known multiferroic BiFeO"3 [19] and is also consistent with the 
strong phonon instability at the R point found in cubic BiCrC>3 [72]. A rock- 
salt-like ordering of Fe 3+ and Cr 3+ cations on the B site lattice of the per- 
ovskite structure was enforced in the calculation; such ordered structures can 
be achieved by layer-by-layer growth in epitaxial films, as has been demon- 
strated experimentally for the case of La2FeCr06 [78]. The calculated sponta- 
neous electric polarization of Bi2FeCrO"6 was about 80 ^C/cm 2 . The calculated 
ground state magnetic structure was the ferrimagnetic equivalent of G-type 
antiferromagnetic ordering, with oppositely oriented magnetic moments of the 
Fe 3+ and Cr 3+ ions resulting in a magnetization of 2 ^/formula unit. Thus, 
ordered Bi2FeCrO"6 was predicted to be a unique magnetoelectric multifer- 
roic with large electric polarization and large magnetization [77]. The mag- 
netic ordering temperature of the newly predicted magnetoelectric multiferroic 
Bi 2 FeCr0 6 , estimated by calculating the nearest and next-nearest neighbor 
coupling constants from first-principles and then applying the mean-field ap- 
proximation, was shown to not exceed 100 K [79]. Since this mean-field value 
can be viewed as an upper bound for the real magnetic ordering tempera- 
ture it is unlikely that Bi2FeCr06 will show multiferroic behavior at room 
temperature. 
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In order to gain more insight into the factors that promote high magnetic or- 
dering temperatures in such ferromagnetic ferroelectrics, a comparative study 
of the series of multiferroics BiFeC>3, Bi 2 FeCr06, and BiCrC>3 was performed 
[79]. The magnetic ordering temperatures of the different compounds were 
again obtained by calculating the nearest and next nearest neighbor exchange 
coupling constants using first principles methods and determining the result- 
ing ordering temperature within the mean-field approximation. It was found 
that the calculated variation of the strength of the nearest neighbor coupling 
constants could be explained based on the different electron configurations of 
the constituent ions using Anderson's theory of superexchange [75]. In addi- 
tion, the influence of the structural distortions on the magnetic coupling was 
determined by comparing the calculated coupling constants for the relaxed 
ferroelectric structures with those obtained for the corresponding undistorted 
cubic structures. In general, the structural distortions resulted in a weakening 
of the antiferromagnetic coupling, and in the case of Bi 2 FeCr0 6 an interesting 
crossover behavior from antiferromagnetic to ferromagnetic nearest neighbor 
coupling could be observed, again in good agreement with the qualitative pre- 
dictions from the theory of superexchange [79]. The insight gained by studying 
the systematic variations of exchange interactions in this series of compounds 
was used to propose a variety of other systems with potentially larger magnetic 
coupling. Since the theory of superexchange is only able to give qualitative pre- 
dictions, however, additional first-principles calculations will be necessary to 
make quantitative predictions. 



5 Summary and Outlook 



From the examples discussed in the present paper it is clear that first princi- 
ples methods are very powerful tools for studying multiferroic materials and 
magnetoelectric coupling. These methods facilitate the accurate quantitative 
calculation of many material properties and can therefore be used to establish 
benchmark values to verify experimentally obtained data. Furthermore, the 
detailed analysis of the calculated results and the possibility to perfectly con- 
trol the system under investigation allows the development of new theoretical 
concepts and the construction of phenomenological models to gain a deeper un- 
derstanding of the basic physics behind multiferroic phenomena. This knowl- 
edge can then be used to design new materials with desirable properties and 
specific values can be calculated explicitly for various candidate materials. 
The examples discussed in the previous sections illustrate how first-principles 
calculations have been used to deliver benchmark data for the magnetiza- 
tion and polarization in BiFeC>3, to understand the mechanism underlying the 
ferroelectricity in YMnOs, to establish the concept of electric- field induced 
magnetization switching in weak ferromagnets, and to design new multiferroic 
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materials with formerly unachievable properties. 

In the future, first-principles methods will undoubtedly continue to make in- 
valuable contributions to the field of magnetoelectric multiferroics. For ex- 
ample, efforts are underway to describe effects of spin-phonon coupling [80], 
and to identify a material that exhibits electric-field switchable weak ferro- 
magnetism [65]. In addition, new classes of multiferroic materials remain to 
be explored, for example the orthorhombic and hexagonal rare earth mangan- 
ites, which show a variety of interesting coupling effects between magnetic and 
ferroelectric properties [8,11]. Nanostructured magnetoelectric composites [81] 
offer a new way of coupling magnetic and ferro- or piezoelectric properties, and 
the multitude of possible materials combinations represents a way to tune the 
properties of the heterostructure to the desired values. Again, first-principles 
calculations offer a powerful tool to study the properties of such heterostruc- 
tures and to predict the characteristics of different materials combinations. We 
therefore believe that the first-principles study of multiferroic materials will 
continue to be very exciting and that these methods will continue to contribute 
significantly to the progress in this field. 
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